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Abstract 

We study a dissipative system of nonlinear and nonlocal equations 
modeling the flow of electrohydrodynamics. The existence, unique- 
ness and regularity of solutions is proven for general L 2 initial data in 
two space dimensions and for small data in data in three space dimen- 
sions. The existence in three dimensions is established by studying 
a linearization of a relative entropy functional. We also establish the 
convergence to the stationary solution with a rate. 
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1 Introduction 



In this paper, we study the following nonlinear system of equations; 



u t + u- Vu + Vp = Au + (1) 

V • u = 0, (2) 

v t + u • Vv = V • (Vv - vVcp) , (3) 

w t + u • Vw = V • (Vw + wV4>) , (4) 

A(p = v — w (5) 



in f2 x (0, oo) for a connected, bounded, open subset of R n with 
smooth boundary dQ. Here u(x, t) is a vector in R n andp(x, t), v(x, t), 
w(x,t) and 4>(x,t) are scalars. Equation (pQ) is the force balance equa- 
tion of a viscous, incompressible fluid with velocity u and incompress- 
ibility condition (|2|). These are coupled with conservation equations 
Q-dH of a binary system of charges with densities v, w and the electric 
potential (ft determined by the Poisson equation ([5]). The force exerted 
by the charged particles on the fluid is A(j)V(j) = V • a(= S™ = i(c r ij)xJ 
where the electric stress a is a rank one tensor plus a pressure; for 
i,j = 1, . . . ,n, 

The electric stress a stems from the balance of kinetic energy with elec- 
trostatic energy via the least action principle, [RLZ07]. For simplicity, 
we have assumed that the fluid density, viscosity, charge mobility and 
dielectric constant are unity 

Solutions for the velocity field equation are determined by the 
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Dirichlet condition 



u(x,t)=0 for (x,t) e dQ x (0,oo). (6) 

Solutions of the equations for the charges are determined by the nat- 
ural (no flux) boundary conditions 

dv d(f) dw d(j) tn . 

— - W -I1 = 0, — +™-^ = 0, ondftx(0,oo), (7) 

ov ov ov ov 

where v is the outward pointing normal to dil. Along with ([6]), equa- 
tions ([7]) are the natural boundary condition of (|3j) and (Q; there is 
no flux of the charges through the boundary. The integral of v and w 
are conserved quantities; 

d r d r , 

v ax = — I w ax = 0. (8) 



dt J n dt Jn 

We assume that 



0< / vodx, wodx<po<oo. (9) 

The constant po is the characteristic charge. The small data results 
in section [6] will in part be formulated in terms of the size of pq. 
Solutions of © are determined by 

(f){x,t)=0 for (x,t) G dn x (0,oo). (10) 

Condition (|10p states that the boundary of the domain is held at a 
fixed potential. The evolution is determined by initial conditions 

u(x, 0) = u (x), V • u (x) = 0, 

(11) 

v(x, 0) = t>o(x) > 0, w(x, 0) = wo(x) > for x G il. 

Since the Navier-Stokes equation is a subsystem of (P)-©, one 
cannot expect better results than for the Navier-Stokes equations. In 



the absence of a fluid, the hydrodynamic system ([I])-© reduces to the 
subsystem 

Vt = V • (Vi> - vVcj)) , 
< Wt = V • (Vu> + roV^) , (12) 
Acp = v — w. 

The equations (|12p are the Debye-Hiickel system, a basic model for 
the diffusion of ions in an electrolyte filling all of R 3 first studied by 
W. Nernst and M. Plank at the end of the nineteenth century, [DH23J. 
The results for the Debye-Hiickel system are complementary to those 
of Navier-Stokes. In |BHN94| . Biler et al. proved the existence of 
global weak solutions in dimensions two and proved the uniqueness 
and regularity of local solutions in all dimensions under appropriate 
assumptions. It is not known whether (|12p in general possesses global 
weak solutions in dimensions greater than two. In view of this chal- 
lenge, one cannot expect solutions of ([I])-© to exist in dimensions 
greater than 2 for general data, c.f. theorems [2] and [5j 

There are several theoreical difficulties associated with the system 
CO)-©. The coefficient Vcfr in the boundary value problem ©-© 
and © is determined by nonlocal information. The right-hand side 
of © and © cannot be formulated in terms of a Frechet derivative 
of a functional (c.f. the Erikson-Leslie theory for liquid crystals and 
the Allen-Cahn and Cahn-Hilliard theory for fluid/interface motion). 
Furthermore, it is not clear how the basic energy law for ©-© (c.f. 
inequality (|18p ) implies the extension property for local solutions in 
dimensions greater than two. Many of the standard techniques for 
parabolic PDE, e.g. the maximum principle and apriori estimates, 
are difficult to apply. 
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Biler et al. presented the first mathematical existence, uniqueness 
and regularity results in [BHN94| . They established the L 2 conver- 
gence of solutions to the stationary solution without a rate. In [BDOOJ, 
this result was improved be establishing an exponential L 1 convergence 
with a rate depending only on Q. Their work relies heavily on the tools 
developed in |AMTU01j and |UAMT00| for the Fokker-Plank equa- 
tion. In [BAMV04J, an exponential L 2 convergence result by means 
of a linearization of the an appropriate energy functional was proved. 

The hydrodynamic setting presented here has been studied in [FG09J, 
[Li09j and [RLW06j . The work of [FG09| establishes several impor- 
tant estimates for the hydrodynamical system when the boundary is 
assumed to be electrically insulated. The work of [Li09] studies the 
interesting zero-dielectric limit of the system on the flat ra-dimensional 
torus. 

1.1 Basic Energy Law 

We develop the basic energy law for electrohydrodynamics. Let us 
consider a classical solution u, v,w of (fT|)- (fTl]) on £1 x (0, T). Assume 
that v and w are positive on Q. Throughout the paper, define 

tpr(s) = s log(s/r) — s + r, r £ (0, oo), s £ [0, oo) 

and define 

^(0 = ^1(0 

The following energy functional will play an important role; 

W = J^(v) +1>(w)dx + i||V0|| 2 + i||u|| 2 2 . (13) 
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The first two terms in this definition are the entropy of the charges 
v and w respectively, while the last two are the electric energy of the 
charges and the kinetic energy of the fluid respectively. 
Differenting W with respect to t gives 

— = (vt,ij/(v)) + K^'H) + (V&, V0) + (u t ,u) 

Here (•, •) denotes the usual L 2 inner product on CI. Integrating by 
parts and using (fTUj) . we see that (S7<pt,S7<p) = — (A0 t ,0) = (wt — vt,(p) 
yielding 

^ = (v t , i//(v) - cf>) + (wt, ij/(w) +4>) + (u*, u) (14) 

The quantities ip'(v) — (j) and ip'(w) + <j> are called the electro-chemical 
potential of v and w respectively. Note that Vt = V ■ (v V log(f e _< ^)) = 
V • (fV(V' / ('u) — </>)) Integrating by parts and using ((7J) 

fa, V> ) - 0) = -(«, |Vlog(t>e^)| 2 ) - (u • Vu, ij/(v) - <P). (15) 

Similarly, by © and ©, 

(wt,il/(w) + <f>) = -(w, |Vlog(u>e^)| 2 ) - (u- Vw,i//(w) + </>). (16) 

Since u is divergence free vector field which vanishes on the boundary 
of Cl, integration by parts gives 

= (u • Vu, u) = (Vp, u) = (u • Vv, ip'(v)) = (u • Vw, ip'(w)). 

Similarly, since (j) is a solution of the Poisson equation ([5]) (u • V(v — 
w),4>) = — (u • V</>, A</>). From ([T]), we have then 

(u t , u) = (Au, u) - (Vp, u) - (u • Vu, u) + (u • V0, A0) 

(17) 

= -|Vu| 2 + (u • V(j),A(j>). 
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Adding (fT5j) . (fT6j) and (fT7|) together the following terms, which are 
interpreted as entropy production due to transport and the kinetic 
energy production due to forcing, 

(u-Vd,^), -(u-Vw,cp), (u-V<p,A<p) 

cancel. We find 

^L = - f u|Vlog(?;e-*)| 2 +u;|Vlog(^e*)| 2 + |Vu| 2 da;<0 (18) 
dt Jo, 

for all < t < T. 

Remark 1. The identity (1181) is the basic energy law for the hydrody- 
namic Debye-Hiickel model. It, along with ([8]), implies that 

H^),^)llLlogL + ||V ( / ) (t)|| 2 2 + ||u(t)|| 2 2 + f ||Vu|| 2 2 ds < 

Jo 

||wo,^o||LlogL + ||V0 O |Il2 + ||u ||l2, t > 0. 

Here 4>q is the solution of the Poisson equation with right hand side 
— ^o- It will be crucial in estabilishing a uniform L 2 estimate when 
dimS7 = 2, c.f. lemmas [5] and [HI 

Remark 2. If we assume that u, v,w and 4> are a classical solution of 
([T|)- (fTH and (j) satifies the boundary condition 

^ = 0, (x,t) € O x (0 5 oo) (19) 
ov 

in place of (|10p or one replaces ©, (J7J), and (|10p with the assumption 

ft = T n (ra-dimensional flat torus), (20) 

then an additional energy law holds; 

d „ „„ 4(p - 1) „„ „ „„ 

T kii) ?„ + — \\Vv,Vw fa < 0, Kp<oo. 



One has equality in the above relation if and only if v(x, t) and w(x,t) 
are equal for a.e. x G Q. Assuming either (fT9|) or (|20j) . a necassary 
condition for a static solution (one where vt = wt = 0, u t = 0) to 
exist is 



The appeal of electrolyte fluids in application is the presence of sharp 
boundary layers in the charges and potential when a static equilibrium 
is reached. Prom the point of view of the analysis and physicality of 
the model, the assumptions (fl"9j) or ([20]) lead to an over simplified 
model. The boundary conditions used in the model in this paper are 
the physical ones but lead to significant difficulties in the analysis. 

1.2 Definitions 

The following function spaces will be used throughout this paper. 
U k (n) = W k ' 2 (n) is the Sobolev space with norm || • || H fc, 



H 1 = dual of H^Q), 
V(n) = Cg°(ft; R") n {v : V • v = 0} , 
H(O) = closure of f in L 2 (ft), 
V(O) = closure of f in H 1 ^), 
V* = dual of V(n), H* = dual of H(J2). 




In this case, the static solution will be 



<p = const, v = w = const. 



In fact, the static solution is the stationary solution 
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The dependence on 0, will be omitted when the context is clear. For 
< t < T < oo, define 



For k a positive integer and < a < 1, f2 C R n and U C R ra x (0, oo) 
open, the Holder spaces 



are defined in chapter 3 of [Lie96| . Note that C fc+a (f2) is defined 
with respect to the Euclidean distance while the space C k+a (U) is 
defined with respect to the parabolic distance. Recall that if k is an 
integer and v G C 2k+a (U), U C R n x (0, oo), then v has 2k uniformly 
Holder continuous derivatives in x and k uniformly Holder continuous 
derivatives in t both with exponent a. 



Throughout, p* will denote the critical Sobolev exponent — Let T > 



0,u G L 2 ((0,T);L«*nV) and Vc/> G L 2 ((0, T); L 9 * ) where ^ + ^ = \. 
We say v G L 2 ((0, T); H 1 ) is a weak solution (c.f. |Lie96] . Chapter 10, 
Section 6) of the equations 

v t + u • Vv = V • (Vv - uV0), x £tl,t> 0, 



Q T = ^x(0,T), Q (li T) = Ox(t ) T). 



C fc+Q (C7), C 



(U) 



1.3 Weak Solutions 



< 



(9u 50 
v(x, 0) = fo(a;) 



x £dQ,t>0 



(21) 



on Qt prodived 




\/uj G H , a.e. < t < T. 
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Setting uj to be the constant function 1, we see that a weak solution 
satisfies the conserved mass equation, (|8j). Note that if v G L p and 
V0, u£L' for a.e. t G (0, T), then vV(j) and uu lie in L 2 for a.e. t and 
hence the trilinear terms are well defined. By virtue of being a weak 
solution, v is absolutely continuous from [0, T] into L 2 . The Sobolev 
space H 1 is separable and by the Lebesgue differntiation theorem, 
v t G ^((O,! 1 );!!- 1 ) exists. 

Let T > and f G L 1 ((0, T); L 2 ). We say u G L 2 ((0,T);V) is a 
(the, incase n = 2) weak solution (c.f. [TemOlj . Chapter 3) of the 
equations 

u t + u • Vu + Vp = Au + f , 
V • u = 0, 

(22) 

u = o, igSD.oo, 

u(x, 0) = uo(x), V • uo = 0, x G $7 
on provied 

(u(t) - u , v) = - / 6(u, u, v) + (Vu, Vv) + (f , v) ds, Vv G V, 

a.e. t G (0,T). 

We have defined 



6(u,v,w) = 




If n < 4, then p* > 4 and &(■,-,•) is well defined on (H 1 ) 3 . If u is 
a weak solution then it is absolutely continuous from [0, T] into L 2 . 
The space V is separable and by the Lebesgue differntiation theorem, 
u t G L 1 ((0,T);V*) exists. 
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If the body force f takes the form A<^V<^>, then we will sometimes 
define 

(f,v) = -(V0®V0,Vv), VuGV. (23) 

The motivation for this definition is easily seen from the identity V • 
(V0<8>V</>) = A(/>V<^+^V|V(/)| 2 and integration by parts. The relevant 
consequences are as follows. Assume <f> solves the Poisson equation 
© and CED. If for 1 < p < oo,n < q v t ,w t G L p ((0, T); H^" 1 ) and 
v,w G L°°((0,T);L*), then by the injections H 1 ^ Hj, L<? wj' 9 
induced by the Poisson equation and the Sobolev embedding W 1,q C 
C a , 

||f t || v * < 2||V0<g> V0 t [| L2 

<2||V0|| L oo||V^|| L 2 (24) 

< c(fi) ||u s wIIl, Hvt, t«t Hn-1 G L p (0,T). 
Let T > and uo G H and Vq,Wq G L 2 . We say (u,v,w,(p) are a 
lueaA; solution of the boundary value problem (fTj)- (|llj) on provided 
u G L 2 ((0,T);V) and v,w G L 2 ((0, T); H 1 ), and for all v G V, cj G 
H 1 , T] G Hj and a.e. < f < T, 

(u(i) -u ,v) = - f 6(u,u,v) + (Vu,Vv) + (A</>V0,v)ds, (25) 
jo 

(v(t) -v ,u) = - (Vv-vV<j)-vu,Vu)), ds, (26) 

-/ (Vw + wX7(f>- wu,Vu), ds, (27) 
■/ o 

-(V<f>(t)-Vr ) ) = (v(t)-w(t), V ) (28) 



A weak solution (u, v, w, <fi) of (fl~1)- (fTTj) is said to be global if it 
is defined for all T G (0, oo). Note that by virtue of the injections 
H fc w H fc+2 n Hq, from the Poisson equation and the Sobolev em- 
beddings H 2 C W 1 ** C C a , if v,w G L 2 ((0,T);H 1 ), then A<£V<£ G 
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L 1 ((0, T); L 2 ). It makes sense to speak of the weak solution u of the 
Navier Stokes equations for then f = Similarly, if n < 4, then 

u e L 2 ((0, T);L q * n V) and it makes sense to speak of a weak solution 
v and w. 

A weak solution is said to be classical if it is possesses enough 
differentiablity to satisfy the equation(s) continuously in the usual 
sense. 



1.4 Stationary Solutions 

One arrives at the stationary equations of (fTj)- (|10p by setting all 
derivatives in t to zero and setting u = 0. They are 

V • (Vv - vVcj)) = 0, 

V • (Vw + wVcf)) = 0, x e O 

Acp = v — w, Vp = A(frV(j), ^29) 
dv d<f> dw d<f> 

v— = — + w— = cj) = 0, x£dn 

ou ou ov ou 



vdx= / vodx, / wdx= / wodx. 
: Jn Jn Jn 

We say v,w G H 1 ,^ G Hj is a weak solution of (f29|) provided the 
equations are satisfied weakly in the usual sense. 

Theorem 1. Let £1 C R n be open, bounded with smooth boundary 
d£l. Let M and N be positive constants. Then there exists a unique 
4> 6 Hq satisfying 

e 4> e —<t> 
= M-p — -7—— - N- r — . (30) 

If k is nonnegative, then cf> G C k (Q) and 

M lim o ||0|| cfe(n) =O. (31) 
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Proof. The unique solution is found by means of the direct method of 
the calculus of variations. For (j) 6 Hq, define 

J[(j>] = J||V0||^2 +Mlog / e*dx + iVlog / dx. (32) 
2 Jn Jn 

Equation (f3"Uj) is the Euler-Lagrange equation of J[-]. Let </>, ^ G L , 

4> ^ ifj, and < A < 1. By Holder's inequality, 

log y e^ 1 -^ dx < log j . ^y ] | 

= Alog / dx + (1 - A) log f dx. 
Jn Jn 

This shows that J is strictly convex. By Jensen's inequality, 

J[<P] > hvm 2 + (M - N) [ cj>dx+ (M + iV) log |fi|. 

Then, by Holder's inequality and the Poincare inequality, J is bounded 
below by some constant depending only on M, N and Q. By the direct 
method of the calculus of variations, J has a unique minimum (j) £ Hq. 
Since J[4>] < oo, it follows that e^, S L 1 and 4> satisfies 



/ 

Jn 



h e * d y Jn e 



By means of the identity, ae t — be 1 = (ab)^ sinhjt — |(^)} a,b > 
0, t E R, this equation is equivalent to 



- / V0 • Vi/> = a[<f>] / sinh(</> - 0[</>])ij> dx, 6 Cg°. (33) 
where a [it] and are defined by the relations 

ui ( MN \* flui X i f N !n e **V \ 



Since sinh(-) is increasing, applying the maximum principle to (|33p we 
find that G L°°. It follows from the standard theory for semilinear 
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elliptic equations (e.g. |Tay97| , Chapter 14), that (ft S C fc (S7) for all 
k > 0. 

If (ft G Hq and J[(ft] < oo, then 

J|(1 _ A) ^ A01 _ ja £ _ 

A 

H (ft 7^ (ft, then the right hand side is negative and (ft is not a critical 
point of J. Hence any solution of (|30p is identically 0. 
Assume without loss of generality that M < N. Define 

C = A(ft. 

One checks that £ satisfies the equation 

AC = (|V0| 2 + /7)C 

where 

= M —p — — — + AT- 



The function r/ is positive. By the maximum principle, £ has no neg- 
ative internal minima nor positive internal maxima. Since 



(dx = M-N < 0, 
we see that £ is nonpositive. Since (ft restricted to d£l is 0, 

(ft(x) < 0, Vx € O. (34) 



Choosing the argument of (|32p to be the constant function 0, we 
see that 

J[(ft] < (M + iV) log 

Arguing as in the beginning of the proof we find that there is a constant 
ci depending only on Q for which 

||V0|| L2 < Cl (AT-M) 2 . 
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Applying Jensen's inequality and the Poincare inequality once more, 

/ eUx> \n\e c ^ M ~^ 2 (35) 
Jn 

for some constant C2 depending only on Q. 
Combining ([Ml) and (l35|) . 



By the maximum principle, <j) converges to uniformly in Q as M, N 
converge to 0. The estimate ([3T]) now follows by bootstrapping the 
elliptic estimates for <fi. □ 

Corollary 1. There exists a unique solution Voq , w^ , (poo of the sta- 
tionary equations // k > 0, then v oo , w^ , (j)^ G C fc (^) a?i<i 

Voo= l v dxe~*°° , Woo = wq dxe^ , 
Jn Jn 

lim ■ [|foo|lc*ccn = lim n ll^oo|lcfe((7) =°- 

po— >0 v ; po— *0 v ' 

Proof. Suppose v, w G H 1 , G Hq is weak solution of (f29|) . For (5 > 0, 
consider the test function 

z/> = log{u + e _ * + £} G H 1 . 



Multiplying the first equation in f|29|) by ip and integrating by parts, 

f f e~^Wv - vV6\ 2 
0= / (Vv-vV(f>) ■ Vipdx= / 1 — -r -^dx. 

This implies that Vv = vV4> for a.e. x £ {y £ Q : v{y) > 0}. If 

i> = \og{(-v)+e-* + 8} GH 1 , 
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then the above reasoning also implies that Vv = vX7<f> for a.e. x G 
{y G 0, : u(y) < 0}. Using the differentiability of Sobolev functions on 
lines (e.g. |EG92j . Section 4.9), we find that 



v(x) = / v dxe for a.e. 

Jn 



Similarly 



w(x) = / wo dxe^ x \ for a.e. x G O. 

It follows that e^, G H 1 . 

Returning to Theorem [H is the unique solution of (|30p with 
M = J* n dx and TV = vq dx. It follows that v and u; are, a fortiori, 
unique. □ 

Remark 3. Let (Too = V(/>oo (g> V^oo — ^|V0oo| 2 /. Observe that Observe 
that 

V • (Too = A^ooV^oo = V(Voo + Woo). (36) 

Equation (136p states that the divergence of the stationary electric 
stress doc is the gradient of a pressure. This is consistent with the 
fourth equation in (|29p . 



1.5 Main Results 

The following existence, uniqueness and regularity theorem is the 
strongest result expected from (fT1)- (fTTT) for general data. 

Theorem 2. If dimQ = 2, and uq G H andvo,Wo G L 2 , then (fB- ffTTI) 
possesses a unique, global weak solution. The solution is classical. In 
particular, ifO<t<T<oo and is any compact subset of Qt, 
then 

ueC 2+t "«nC a (Q (i]T) ), v,w G c 2+a (Q (tiT) ). 



16 



Thus, in two space dimensions, (fTj)- ([TTj) is solvable and enjoys usual 
regularization property found in equations of parabolic type. In di- 
mensions three and four, the existence of a global weak solution can 
be proved using the techniques in the proof of theorem [2] assuming a 
uniform in time L 2 apriori estimates for v and w. 

The next theorems concern the long term behavior of weak solu- 
tions. In order to quantify the convergence, define for 1 < p < oo, 

S p (t) = 

I , / m9 \v(t) — t>ool P \w(t) — U?ool P i __ ( . < „ , .9 , 

/ |u(t)| 2 + 1 W °°' + ' ^J™ 1 +|V0(t)- V^oo| 2 dx 

The following theorem is modeled after [BDOOj . Note that, in contrast 
to theorem [21 these solutions may not be defined globally if n > 3. 

Theorem 3. Let O C R n , be bounded and uniformly convex, n > 2, 
and {u,v,w,4>) be a global weak solution of (fT])- (jlip . Then there is 
Ai > depending on f2 and e\ < oo depending only on the initial data 
so that for all t > 

£\(t) < ei e~ Xlt . (37) 

It is difficult to extend the technique of [BDOOJ used in the proof of 
theorem [3] to Q with general geometry. However, if one assumes that 
the initial data is close to the stationary solution and the stationary 
solution is small then the convergence with a rate is recovered. This 
result is a kind of linearization of the argument used in the proof of 
theorem [3j 

Theorem 4. Let SI C R n , n = 2,3, u G H,vo,wo £ L 2 . There are 
positive constants 

P2 = P2(ty, A 2 = A 2 (T2), e 2 = e 2 (£l) 
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such that if 

<%(0) < e 2 , po < P2 
then (fT])- (|lip possesses a global, weak solution and 

Mt) < e 2 e- X ' 2t . (38) 
If a global weak solution satisfies 

sup \\v, 1L>|| L 2 < oo, 

te(o,co) 

f/ien i/iere is to > so that ^2(^0) < e 2- 

Thus, in two space dimensions, theorem [J] implies that the solution 
from theorem[2]tends to the stationary solution since the weak solution 
is eventually close to the stationary solution. 

In three dimensions, a global existence, uniqueness and regular- 
ity result is proved under a small data assumption. However, it not 
sufficient to assume that the initial data is small and close to the sta- 
tionary solution. One must also assume that the stationary solution 
is also small. 

Theorem 5. Let dimQ = 3, and uq £ H and t>o, wo £ L 2 . There exist 
constants 

p 3 = p 3 (Q), e 3 = e 3 (Q), 5 3 = 5 3 (Cl) 

such that if 

po < p 3 , £2(0) < e 3 , ||u ,wo,wo||h2 < 03, 

then (fTj)- (lllj) possesses a unique, global classical solution. In particu- 
lar, if < t < T and & is any compact subset of Qt, then 

u e C 2+a (3?) n C a (Q (m ), v,we C 2+a (Q (tiT) ). 
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The remainder of the paper is organized as follows. In section O 
some sufficient conditions are developed for concluding L°° bounds on 
weak solutions are developed. These are later used proofs of theorems 
[2] and EJ As noted in the introduction, this is a necassary step for the 
regularity programme due to the lack of a gradient descent structure 
for the conservation equations and the lack of a maximum principle for 
Q and (HD . In section [U local weak solutions are constructed and the 
extension property is developed in two space dimensions. The section 
is concluded with the proof of theorem [2j In section [5l we give a proof 
of theorem modeled after the result of |BD00j . Finally, in section 
we give several preparatory lemma and provide the proofs theorems U] 
and[5j 

The letter C will denote a constant which may change from line to 
line within a proof. The letters c\ , oi , . . . , will denote constants that 
are fixed throughout the paper. 
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3 Preliminaries 



In this section we are assuming Q, is an open subset of R 2 or R 3 and 
that <9S7 is smooth. The following preliminary results will later be used 
to infer uniform L°° bounds on solutions v, w, u. The C 2+a regularity 
of solutions will then follow from classical results for linear second 
order PDE of parabolic type. 

In the case of L°° bounds on weak solutions v and w of ([3]) and 
([4j), the usual techniques for nondivergence form semilinear parabolic 
equations do not apply, mainly due to the boundary conditions ([7]). 
Instead, we rely on a Moser type iteration argument. The essential 
part of the argument is that due to the divergence free condition, no 
regularity on the velocity u need be assumed. 

First a general 

Lemma 1. For t G [0, T] and p £ [l,oo), let y(t,p) be positive and 
continuous and satisfy the differential inequality 

^(t,p)+yP(t,mp) < w(t)p k y a ^(t,p) 

where m > 1, k > 0, w(t) > is measurable with w(r) dr < j\t — s\@ 
for some constants 7, (3 > and a(p) G (0, 1]. IfO < e < t and 1 < po, 
then there is C3 = c^(w, k,po, e) < 00 for which 

lim y(t,p) < c 3 • y(t - e,p ). 

pjoo 

Proof. Without loss of generality, we may assume that a(p) = 1 for 
every p. For otherwise, we may replace y(p,t) by max{l, y(p, t)}. Let 

< s < t < T. Using Gronwall's inequality, 

"t / ft 



y p (t,p) + J y p (r,mp) < exp (^p k j w(r)dr^j y p (s, 



p) (39) 
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We will take advantage of the various powers in this inequality to infer 
some bounds. 

Let < e < t and a > 2. Define 

U = t- ea~\ Si = t i+ i - ti, pi = mVo, i = 1,2,... 

From ([39]) with p = Pi+i,t = and s £ [ij, tj + i] we have (recall 
mpi = p i+1 ) 

( Pi+i f u+i \ 

2/ Pl (tj + i,p i+ i) < exp — — / w(r)dr y Pl {s,mpi). 

V m Ju J 

Integrating this expression with respect to r = s G 

( Pi+i f ti+1 \ f U+1 

5iy Pl (ti +1 ,p i+ i) < exp I -!— J w(r)dr\ J y Pl (r,mpi) dr. 

Using ([39]) once more with s = ti,t = ti + % and p = pi to bound the 
second integral on the right hand side, 

6fy*(t i+1 ,p i+1 ) < exp (^(^±L + p^j w[r)d^j y p *(t uPl ). 

This in turn implies 

y(t i+1 , Pi+1 ) < exp (( ^ W ( r ) dr j 

Note the integrand in the argument of the exponential is bounded 
above by 

Clearly 2 ecr~* < <5i < e<7~\ Choosing <r = max{2,m z 3 }, 

R < Cim~ kl 

for some C\ = Ci(e, 7, (3,po,m, k). Similarly, 

- — -L 
Si w < C 2 {e,P,p ,m)^ . 
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It follows that for all i = 1,2..., y(t i+1 ,p i+1 ) < e ClTn kt y(U,pi) 
and so by recursion 

y(U+i,Pi+i) < $iy{t ,p ) 

i 

where $j = n* =1 C^™* exp(Cim" fe *). This product converges and we 
set the limit to be C3. The conclusion now follows by varying po and e 
in a sufficiently small subset of [0, T] x [1, 2] and using the continuity 
of y. □ 

Proposition 1. Let S > 0, dimfi = 2,3, 

vGL 2 ((0,5);V), 

and vo,wo G L 2 . 

L Then there is < Tq = To(S7, ||t>o, ^o||l 2 ) < S so that the 
problem 

v t + v • Vt> = V • (Vw - v\7(f>) , 
u> t + v • Vu> = V • (Vu> + wV<f>) , 
= v — iw, 

£-1^ = 0, ^ + 4^ = 0, = 0, O n^x(0,oc), 
ai/ az/ Of £71/ 

t> (x, 0) = vo(x), w(x,Q)=wo(x) x£fl. 

/ias a unique weak solution on Qr - Moreover, if p > 2, i/iere is a 
constant C4 = 04(0, sup tg ( ,T ) II u )HIl 2 >t?) so 

sup ||v,u>||lp < ||uo) ^o||Lpe C4T . (40) 

*6[0,T] 

2. Ifp > 2 and < s < t < T , then there is 

C 5 = C 5 (n, ||w,w|| L 2((0,T);Hl) ) P ) i - S) < OO 
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for which the weak solution v, w satisfies 



?;(*), < c 5 • \\v(s),w(s)\\ LP . 



(41) 



3. If, additionally, 



v , w G C 2+a (Q), < v (x),w Q (x) Vi6Q, 




(42) 



on dQ, 



and v G C a (QT ), i^en 



v,w E 



(QT ) 



and f and u; are positive on Q y . 

Proof. (Part 1) The existence and uniqueness of a weak solution is 
established by making slight modifications to the proof of |BHN94j , 
theorem 1 to account for the term v • V. We omit the details here. Let 
< To = Tq(U, ||vo||l 2 ) II^oIIl 2 ) so that the weak solution is defined 
on Qt and let 



Multiply the v -equation and the w-equation by pv p and pw p ~ 
respectively and integrate over f2. Since pVvv p ~ l = X7v p andpVwu^ -1 = 
Vw p , using V • v = and v £ Hq = we have 



Mo = SUp H^HIl 2 ; M\ 

te(0,T o ) 





da- 




Integrating by parts gives 




4(p - 1) 



Vv 2 , Vw^ 




V 



(43) 
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for a.e. t G (0, To). Since 4> solves the Poisson equation, by elliptic 
regularity and the Sobolev embedding H 1 C L 6 , we have 

||V0|| L 6 <C(n,Mo), a.e. te[0,T ) 

For e > we have from the Sobolev inequality 

p p -i — p p 

\\v2 ,w2 [| L 3 < C(£l)(e~ \\v,w\\£ p + e||Vu2 , Vwa || L 2). 

Applying these inequalities, the estimate 

A<2(p- 1)|| V<£|| L 6 Hut, || L 3||VV2, Vw2|| L2 

< C(Q,M ,p)(e~ 1 \\v,w\\l P +e||Vw5 ) Vw*|| L a)||Vu5,Vu;5||i 1 a 
follows. Choosing e < (2C(0, M ,p)) _1 , we find from (03]) with p > 2 

^lk™llL> < C(n,M ,p)\\v,w\\l P . 

Gronwall's inequality now implies (|40p . 
(Part 2) Returning to (|43p . the inequality 

^IK^IIL + ||V^/ 2 , W/ 2 ||2 2 < 4^1^011^11^,^11^ 

holds for q > 2. By elliptic regularity and the Sobolev embedding 

W 1 ' 4 C L°°, we have ||V0|| L °° < C[Q)\\v, w\\ H i G L 2 (0,T ). With 

this constant, define w(t) = 4C(Sl)||-u(t), ^(i)|| H i and note that 

* i 
w(r) dr < j\t — s| 2 

where 7 = 7(0, Mi). Furthermore, by the embedding H 1 C L , 
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,v, w \\i 2q = 11^,^111,4 < c(n)(\\vv^ 2 ,vw^ 2 \\ 2 ^ + \\v,M\h)- 

Combining this with the previous observations, one has 

^IKHIL + IkHlL^ ^ q 2 w{t)\\v,w\\l q . 
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We may thus apply lemma Q] with y(t, q) = \\v, w||l8, ol = l,j3 = \,p = 
k = m = 2 to find (|4T|) . 

(Part 3) If fo,^o G C 2+a (f2), then certainly ^0,^0 G L p for p > n. 
Fix n < p < p* . Then elliptic regularity and the Sobolev embedding 
W llP C L°° and (JH imply there is C = C(0, M , T ,p) < 00 so that 

||V0|| LO o (QTo) < C, Vt€[0,T ). 

If v G C q (Qt| ) ), then w and w solve divergence form equations with 
oblique boundary conditions with bounded, measurable coefficients. 
By [Lie96| . theorem 6.41 and 6.44, v,w G C Q (Qr ). The relation 
A<j) = v - w then implies V0 G C a (Q To ) as well. The C 2+a (Q To ) 
regularity of v , w is then guaranteed by the smoothness and consis- 
tency assumption (|42|) and [Lie96], theorem 5.18. 

One readily checks that || min{0, v, w}\\i J 2 = for all t G [0,To]. If 
v is not bounded below by a positive constant on Qt > then v(x, t) = 
for some (x,t) G Qt U (dQ x (0,To)). If G Qt , then the strong 
maximum principle ( |Lie96j . theorem 2.9) implies v = 0, contradicting 
(jSJ). If (x, t) G <9f2 x (0, T ), then by the parabolic Hopf lemma ( [Lie96] . 
theorem 2.6) implies 

dv dv dd> 

0>—(x,t) = —{x, t) - v(x, t)^-(x, t) = 0, 
dv dv dv 

again a contradiction. Analogous arguments apply to w. □ 

Solutions of the velocity equations will be constructed from a spe- 
cial basis. To describe the method, we choose a particular orthonormal 
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basis of H satisfying 

AQ + Vpi = -XQ 

v • Ci = o 

(i(x) =0 for x G aT2. 

Here Q £ Cq 3 (Q) and < Ai < A2 < . . . are eigen-pairs of the 
Stokes operator on £1. For i,j = 1, 2, ... , the functions Q satisfy the 
orthogonality relations (c.f. section 2.6, [TemOlj ) 

(Mi ■ Cj) = -Aj(Ci,G) = -Ai<%. (44) 

Following proposition[TJ we are also interested in L°° bounds of the 
velocity u. The following general result for the Navier-Stokes equations 
will be useful in the regularity proof in two space dimesions. 

Lemma 2 (Ladyzhenskaya's Inequality, c.f. [TemOl], Ch. 3). Let O 
be an open subset 0/R 2 or R 3 . Ifu £ V then 

1 1 

|| u IIl 4 — II u IIl2 ||Vu||^ 2 ) dimfi = 2, 



1 3 



II u IIl 4 — II u IIl2 II^ u IIl 2 ' dimfi = 3. 
Proposition 2. Suppose that T > 0,dim Q = 2,uq £ L 2 , 

f £ L 2 ((0,r);L 2 ), t 2 i t £ L 2 ((0,r);V*), 
and u is the unique Leray-Hopf solution of 

— + u • Vu + Vp = Au + f 

at 

div u = 

on Q-p. T/ien there is a constant c§ = c$(p,, f, ft, Uq) < 00 so i/iat 

II l|2 / C 6 ||V7 ||2 / C 6 / ||V7 l|2 j ^ c 6 

||u t || L2 < ^, ||Vu|| L2 < -, y ||Vu t || L2 dr < -j. 
for allO < s <t< T. 
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Remark 4. In proposition [2] as in theorem [21 we are merely assuming 
uo £ L 2 . Also, ft may be singular at the origin. We conclude that u 
is C a away from the set £1 x {t = 0}. 

Proof. For the moment, assume ft 6 C°([0, T]; L 2 ). Let m be a posi- 
tive integer and v(x,t) = Y^=x d( x ) v i (t) an d v% solve m-dimensional 
system of ordinary differential equations 

m 

Vi + XiVi + B ijk VjV k = Fi(t). (45) 

We have defined 

5ijfc = 6(Cj»Cfc,Ci) F i (t) = (f,C < ). 

It is well known (c.f. [TemOlJ, Chapter 3) that for some subsequence 
of m, the v converge to u in the strong topology L 2 (Q^) and there is 
Co depending only on uo and f for which 

II V IIl-((0,T);L2) + II Vv IIl2((0,T);L2) ^ C 0> Vm > °" ( 46 ) 

Differentiate (|45f) with respect to t and multiply the resulting system 
component- wise by Vi. We find that 

(v tt , Vt) + 6(v t , v, v t ) = -(Av t ,v ( ) + (ft, v t ) 

Multiply this equation by t 2 and define w = tVf One easily checks 
that 

^^■II w IIl 2 + II Vw IIl 2 = K w > u > w ) + t ll v *llL 2 +*(ft,w). 

On the other hand, if we multiply (|47p by -Uj componentwise, we find 
||v t || 2 2 = (Av,vt) - 6(v,v,v t ) + (f, v t ). 
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Inserting this expression into the previous equation and noticing that 
t(Av,v t ) = -t(Vv, Vv t ) = -ti|(Vv, Vv) + |(Vv, Vv), gives 

~ [Mh + tllVvl^) + ||Vw||2 2 = A (t) + B(t) + C(t) 

where 

A(t) = |||Vv||i a , 
B(t) = Bi(t) + B 2 (i) = 6(w, v, w) - 6(v, v, w), 
C(i) =t(f t ,w) + (f,w). 

From (I46p . is clear that A(i) is integrable with integral bounded 
above by |co along any measurable subset of [0, T]. Similary, one has 

G(t) < 2i 2 ||f t ||^ + C(n)||f || 2 2 + ~||Vw|£ 2 

with the two left hand terms on the right hand side of the ineqaulity 
also integrable. By Holder's inequality, 

Bi(t) = — / w • Vw • v dx < ||w|| L 4||v||l4||Vw|| L 2. 

•/!! 

i I I 

By Ladyzhenskaya's inequality, ||v|| L 4 < 22 ||v||^ 2 ||Vv||^ 2 which lies 

1 I s 

inL 4 (0,T). Similarly, ||w|| ||w|| L 2 < 2a ||w||£ 2 ||Vw||£ 2 . Taking these 
elements into consideration, we find 

Bi(t) < 8||w|| 2 2 ||v||£ 4 + 7l|Vw||^ 2 . 

Similar ideas imply 

B 2 (t) < 2||v||£ 4 + i||Vw||2 2 . 
In total, we find that 

jF(t) + G(t)<B(t) + E(t)F(t) 
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where 



D(t) = 2t 2 \\{ t \\ 2 / , + 2||f|| 2 / , +4||v||^4 G ^(O.T), 
E(i) = 16||v||£ 4 G 1^(0, T), 

F (0 = II w IIl2 + *II Vv IIl2, 
G(t) = ^||Vw||2 2 . 

We have assumed that fj is continuous from [0, T] into L 2 . Therefore, v 
and vt remain bounded in H 1 and F(i) = [|w||? 2 + ^l|Vv|| 2 j2 = when 
t = 0. By Gronwall's inequality, there is a constant cq independent of 
m and the continuity of fj for which 

F(t) + / G(r) cir < c 6 , V0 < s < t < T. 

J s 

The proposition now follows by first approximating f by functions 
continuously differentiable in t, and then letting m diverge along the 
given subsequence; for a.e. t G (0,T), Fatou's lemma implies 

t 2 \\u t \\l 2 + t||Vu|| 2 2 + t 2 f ||Vu t ||L2 dr 

J s 

ft 



< Hm m (F(t) + G{r) dr^j < c 6 . 



□ 



4 The Local Existence and Extension 
Property in 2D. 

We will construct weak solutions of (fT1)- (lll|) as the limit of modified 
Galerkin approximation. We look for solutions of ([I])-© with the 
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form 

m 

u m (x,t) = y~]uj(t)Ci(x). 

i=l 

The orthogonality of the Q lead us to the following approximation 
problem. For 

Fi(t) = / A(f) m \7(t) m ■ Ci dx, 

consider a solution of 

m 

Ui + Xm + B ijk UjU k = Fi(t) (47) 
j,k=i 

— ^ + u m • Vf m = V • (Vv m - v m V(j) m ) , (48) 
at 

— ^- + u m • Vw m = V • (Vw m + w m V(p m ) , (49) 

= v m - w m , (50) 

Vi(0) = u -&dx, v m (0, •) = v o (0) t« m (0, •) = wo(-)- (51) 
./n 

We will prove 



Theorem 6. Suppose dim f2 = 2, uo £ H and (|42l) ZioZd. For any m > 
0, T G (0, do), t/ie problem (147!) -(|5T j) /ios a unique, classical solution 
with D^D l x u m G C q (Qt) /or A; = 1, 2 and Z = 1, 2, . . . , and v m , w m G 
C 2+a (Q T ). 

Theorem [6] will be a consequence of the following lemmas. 



Lemma 3. Suppose dim $7 = 2,3, uo £ H and assume ([42 p . There 
is To = To(f2, ||uo||l2, H^oIIl 2 ) II^oIIl 2 ) > snc/i t/iai iZie problem (|4"7|) - 
(f5"Tj) Zias a unique, classical solution with D^D x u m G C a (QT ) /or 
A; = 1, 2 and Z = 1, 2, . . . , and u m , u; m G C 2+a (Qr ). 

Proo/. Let s > and let {vj™! G W 1 - 2 ^, s); R m ). Define v(x,t) = 
YT=id{x)vi{t). The embedding W 1 ' 2 ^,*) C Cl(0,s) implies 
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v G Lip(Q s ) (with respect to the parabolic distance.) From proposi- 
tion [H there exists To = To(f2, H^oIIl 2 ) II^oIIl 2 ) an d a solution v,w £ 
C 2+a (Q-r ) of the problem in proposition[TJ Let Mq = sup tg ( 0j r ) ll v o> ^o||l 2 - 
Define Fi(t) = J n A0V0 • Ci dx. Then, by elliptic regularity and ([^0]) 

m m 

sup Vif(t)< sup ||AW|| L i < C(m,n,M ). 

te[o,T ] - =1 te[o,T ] i=l 

Substituting A0V0 into (|47j) . this ordinary differential equation has a 
W 1 ' oo (0,T 2 ) solution {vi}^ v Infact, by flU]) and dSJ), we have 
S^Li <? 4 2 (0) < l u oli- Therefore, the satisfy the estimate 

m 

5^^ 2 (i) < ||u || 2 2 + tC(m,n,M Q ), yt G [0,T ]. 
i=l 

Fix s = To and M = ||uo||? 2 +2oC(m, Q, Mo). This construction maps 
W 1 ' 2 ((0,T );R m ) into 

S M = Jwr=ieW 1 ' 2 (0,To;R m ): sup f^v^t) 2 < M) . 
I i=i J 

The mapping is clearly continuous and takes Sm 1 a compact and con- 
vex subset of W 1 ' 2 ((0,T );R m ), into itself. A fixed point {vi}^ is 
thus guaranteed by the Schauder fixed point theorem ([Lie96], theo- 
rem 8.1.) Let {ui}^L 1 be such a fixed point. Let v m ,w m , 4> m be the 
solution associated with v = Y^IL\ d u i- from this construction. Define 
u m = v = }^ i=1 Qui. 

Note that the L p norms define the same topologies for velocities 
u of the form Y^iLi& u i- Suppose u,v,w and u, v,w are two such 
solutions. Define 

U = U — U, V = V — V, W = W — W, (f> = (J) — (f). 
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Note that v is a solution to the problem 

v t + u • X7v + u • Vv = V • (Vu - uV0 - vV4>) 

dv dtp _d<f> , _ 

■7; W7; t> -7— = 0, on aS2 x (0, 00), v(x,0) = 0, x G SZ. 

dv dv ov 

Define <5(t) = ||u, u, -tD || ^ 2 - Multiplying the first equation by v and 
integrating by parts we find 

j t l\Mh + \\Vv\\h = m),Vv). (52) 

for O(i) in terms of the two solutions. In each of the bilinear terms of 
0, there will be one cross term involving the difference of the solutions. 
Then 

0(i) < C*(t) 

for some constant C depending only on Q, m and the initial data. 
We omit the details as they follow in a straightforward way using the 
equivalence of the L p norms on the velocities, (|40p and the regularity 
property of the Poisson equation. An analogous estimate applies to 
w. 

It is also clear from (I47p . that there is C depending only on fl,m 
and the initial data for which 

^||u|£ 2 <C||u||£ a + 11511^ (53) 

where S = A0V0 — A^V^, and another such C such that 

mi < cm- 

Adding ([52]) to the corresponding inequality for w and to ([53]) . we 
have shown that there is C depending only on Q, the initial data and 
m for which 

j t m < cm 
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Since $(0) = 0, we infer 




uniqueness is now established. 

The C 2+a (Q To ) regularity of v m and w m now follows from propo- 
sition [1] and that fact that u m G Lip(Qr ). The regularity of u m 
follows from the easy observation that the right hand side of ()47|) is 
Holder continuously differentiable in t and that D l x Q S C Q (Qr ) for 



Lemma 4. Let ft be an open subset o/R n ,n = 2,3. If u G H 2 (ft) n 
Hj(fi), tfien 



Proof. We prove the 2-dimensional case first. Without loss of gener- 
ality we may assume u S C^°(R 2 ). Let v = u x . Then v has compact 
support. Clearly 



all I > 0. 



□ 




dim ft = 2, 



dim ft = 3. 




Thus 




Similarly, 




Then 




v(x,y)dxdy< / f(x)dx g{y)dy. 
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However, by Holder's inequality and the relation u = v x , 

/ f(x)dx<3 \v\ 2 \Vv\ dxds < 3||u x |£ 3 ||V 2 u|| L 3. 

Bounding the integral of g in the same way, and using the interpolation 
||Vn||^ 3 < || Vn|| L 2 ||Vu|| L 6, we find 

||«x|Il6 < 9 ll^|lL3||V 2 ti||L3 < 9||u x ||l2||«s|Il6||V 2 u||l3 



or 



\u x \\\fi < 3||Vu|| L 2||V 2 ii|| L 3. 



Arguing similarly with u y , we have, by the triangular inequality, 
||V«|| L 6 < (K||£ 8 + ^ ^6||Vn||| 2 ||V 2 n||£ 3 . 



This proves the first part of the lemma. 

Now we consider the 3-dimensional case. For z £ R, 

I|V<,-,2)Hl3(r2) = JJ^\Vu\ 3 (x,y,z)dxdy 

= 3 / (| Vu| Vu • Vu z )(x, y, s) ds dx dy 

J Jr 2 J-oo 

<3 JJJ \Vu\ 2 \V 2 u\dxdydz < 3||Vu|| 2 3 ||V 2 u|| La . 
By what was shown in the 2-dimensional case, 
JJJ \Vu\ 6 dxdydz < 

72 / \\VU(;-,Z)\\ 
J R 

Consequently, 

||Vit||£ 6 < 72 • 3||Vu|| 2 j3 ||V 2 u|| L 3X 

J ||Vn(-, -, ^)|| L 3(R2)||V 2 n(-, -, «)||l3(r2) dz 

< 6 3 ||Vn||3 3 ||V 2 n||^3. 
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Interpolating L 3 C L 2 n L 6 again, 

||Vu||£ a < 6 3 ||V U ||| 2 ||Vn||^ 6 ||V\|| 3 3. 

Dividing and taking the appropriate power now gives the second part 
of the lemma. □ 

Lemma 5. Let £1 be a smooth bounded subset o/R n and p = and 
q = \p. If there is a constant F so that 



|LlogL < F 

and e > 0, then there is c-j = cj(Vi, n, e, F) for which 

ll/IIL<^llv/ll^ + (i + ll/ll^)c 7 . 

Proof. Without loss of generality, assume / is nonnegative. Let E : 
H 1 (r2) — > H 1 (R n ) be a continuous extension with SuppE 1 / C B(i?) 
and ||.E/||LiogL < F where F = F(Sl,n,F) < oo. Denote / = Ef. For 
M > 0, 

f f q dx< [ f q dx < \B(R)\M q + [ {f-M)\dx 
By the Gagliardo-Nirenberg-Sobolev and Holder inequalities we have 
/ {f-M)\dx<C{n)(f {Dflf^dx) 

JR" \J{f>M} J 



<C(n)\\Df\\l HRn) [J f I** 



Note that 



/ fdx < — !— / f log fdx < 
J{f>M} logM J {f>M } logM 
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Combining the above and applying the continuity of E, 
[ (/ - M)%dx <C(n,n 



F 

1 a„ s nin „\ I r n r\\P 



Cetainly, by the Poincare inequality ||/|| H i < C(fi)(||D/|| L 2 + ||/|| L i). 
Combining this with the results from above 

V 

ll/IIL < c(n,n) (^£j^J 2 \\Df\\i 2 + c(n,n,M)(i + 

Choosing M sufficiently large, the lemma follows. □ 

Lemma 6 QGT01], corollary 9.10.). If Q is a smooth, bounded, open 
subset of~R n , u G W 2 ' p n Hq and 1 < p < oo, then 

||V 2 u|| L p < C(n,p)\\Au\\ hP . 

Lemma 7. // dim Q = 2 and there are constants 01,02,03 so that 

||V0]|l2 < 01, IK^HLiogL < a 2 , ||v,u/|| L i < a 3 , 

E H 2 n Hq solves Acfi = v — w and e,5 > 0, i/ien i/iere is eg = 
cg(^, 6,5,01,02,03) urai/i 



V</> • (fVw — wVw) dx 



< e\\Vv, Vw||l2 — w|| L 2 + eg. 



Proof. By lemma [H with dim $7 = 2, and lemma [6l there is C\ = 
Ci(f2,ai) for which 

||V0|| L6 < 65||V0||J 2 ||V 2 0||£ 3 < CxIKHI^. 

By lemma O with dim £1 = n = 2, there is C2 = 6*2(02, 03, C7) so that 

3 _ 
ty||^ 3 < ei||Vv, Vw|| L 2 +C2. 
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If 62 > 0, then by Nash's inequality, there is C3 = 6*3(0, e 2 , 0,3) so that 
\\vM\l* < £2\\Vv,X7w\\l 2 +C 3 . 

Combining these, we estimate the trilinear term. By Holder's inequal- 
ity, 



/ • (vVv — wVw) dx 
Jn 



< Ci(ei||V«, Vw|| L 2 + C 2 )\\Vv,Vw\\v 

< (C iei + e 3 )||Vt;,Vu;||2 2 + e3 - 1 (2C 1 C 2 ) 2 . 

Choose e% and £3 so that C\t\ + 63 = e/2. Also, choose e 2 = 5~ 1 e/2. 
Then 

^\\Vv, Vw\\ 2 L2 < e\\Vv, Vw\\ 2 h2 - e 2 d\\Vv,Vw\\l 2 

< e\\Vv,Vw\\ 2 L2 - 8\\v,w\\h + SC 3 . 

The lemma follows with c 8 = 5C 3 + 1 (2C 1 C 2 ) 2 . □ 

Lemma 8. Let dim = 2 and let u m , w m , n; m 6e the solution obtained 
from lemmalM For r > 0, Zei tp(r) = r logr — r + 1 and define 

f 11 

V!/ = / i>{v m ) +4>(w m )dx+ — |[V^ m |lx,2 + -H u m|lL2- 

Jn 11 
then there is eg = 09(0,^(0), H^oIIl 1 ) II^oIIl 1 ) such that if 

\\vo,wo\\l 2 < c 9 , 

then 

SUp \\v m , W m \\1 2 < Cg. 

te(o,T ) 
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Proof. From lemma [TJ v m and w m are positive on Qt - Mimicking 
the calculation of the basic energy law (|18p , differentiation of W with 
respect to t gives 

^T = - ! v m \V\og{v m e~^ m )\ 2 + w m \V\og(w m e^ m )\ 2 + \Vu m \ 2 dx 
Jn 

< 0. 

(54) 

for all < t < To. Here we used the identity 

((u m ) t ,U m ) = -||Vu m ||L2 + (Um • V0 m , A(j) m ). 
Recall that the intregral of v m and w m are conserved quantities. 

||^m||r,i = II^oIIli, ll^mllLi = II^oIIli- ( 55 ) 
Then (i54|) implies 

||^ m ,^m||LlogL + l|V^ m |li2 < W(0). (56) 

Define u>(t) = \\v m ,w m \\^j2 and ((t) = \\Vv m , Vw m ||^2- Let e = 8 = |. 
If eg = 2c$, then by lemma (|43p with p = 2, 

^-uj + u + (<c 9 (57) 

This inequality implies u;(t) G [0, max{cg, u;(0)}] for all i G [0, To]. □ 

Proof of theorem [6]. Define 

cio = max{c 9 , W(0), ||u , v , ^o|Il2}- 

The energy decay ([M} implies that ||u m ||£ 2 (i) < W(t) < W(0) < c w 
for all t G [0, To]. The conclusion of lemma [8] implies that \\v m , WmW^it) < 
cio for all t G [0, To]. The solution obtained from lemma [3] enjoys 

||u m ,v m ,«; m ||L2 < 2ci , V£G[0,T ]. (58) 
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Returning to lemma O the solution may be extended to an interval 
[0, Tq + 5] where 5 depends only on c\q. The conclusion of theorem [6] 
now follows from repeated application of this extension property; the 
regularity of u m , v m and w m follows as in lemma [3] because (f22j) is 
satisfied by u = u m (T ),v = v m (T ),w =w m (T ). □ 

Lemma 9. Let dim = 2, T G (0, oo) and u m , v m ,w m be the solution 
obtained from theorem [3 Then there is c\\ < oo depending only on 
the L 2 norms of uq,vo,wq and W(0) but independent of m so that 

f T 

sup \\u m ,v m ,w m \\ L 2 + \\Vu m , Vv m , Vu; m ||L2 ds 
te(o,T) Jo 



+ 



f 




2 


dv m 


dw m 


2 


/o 


dt 


+ 

V* 


dt 


dt 


H" 1 



(59) 



ds < c\\. 



Proof. One can construct a bound of the first two terms in (|59p by 
simply integrating (|54p . (|58p and (|57p. In particular, the L 2 norms of 
are bounded in t uniformly in m. From (|24|) . with p = 2,k = 0, 
||fm||v* is then also bounded in t uniformly in m. It is then standard 
(c.f. [TemOlj . theorem 3.1) to check that 



du r 



dt 



£ L 2 (0,T; V* 



with norm independent of m. 

It is straightforward to check, using lemma [2] and the Sobolev 
embedding H 1 C L 6 C L 3 , that 
dv m 



dt 



= sup (Vf r , 
h- 1 !/eHM|y||i=i 



«mV(/) m - V m U m , Vv) 



< (IbmllH 1 + ll u m||L 4 ll' u m||L 4 + \\ v m ||l 3 II V</)|| L 6 ) . 

Using the Sobolev inequality, we can easily bound the right hand side 
by 

C(lbm||Hi + IbmllHi) 
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for some constant depending only on f2 and c\q. Applying similar 
reasoning to w m , we see then that 

dv m dw m 



dt dt 

with norm independent of m. 



e L 2 (0,T;H -1 ) 



□ 



Lemma 10. Let dimO, = 2, T £ (0, oo) and u m ,v m ,Wm be the solu- 
tion obtained from theorem 0. Let 

f TO = A<f> m X7(p m . 

Then there is c\i < oo depending only on the L 2 norms of vlq,vq,wq 
and W(0) but independent of m so that 



sup ||f m || v * + 

te(o,T) 



f 


dim 


10 


dt 



ds < C\2- 



Proof. (|60l) is an easy consequence of flM]) and (f59l) . 



(60) 
□ 



4.1 Proof of Theorem [2] 

Let uo € H,i>o,u>o G L 2 . Let {vq,Wq} be a sequence of functions 
satisfying (|42j) with t>o> — * v 0' w o i n L 2 as e j 0. Let T G (0, oo) and 
to be a positive integer. Let u m ,v m , w m be the solution obtained from 
theorem [6] with initial data vlq,Vq,Wq. ^ ne bound (|59|) is independent 
of m and e. Applying Lion's compactness theorem ([TcmOlJ, theorem 
2.3), for some subsequence of {u m ,v m ,w m }mioo,tio, 

u e m ^u weakly in L 2 ((0,T); V), weak-* in L°°((0, T); H), 
and strongly in L 2 ((0,T);H), 
v m Mm -+v,w weakly in L 2 ((0, T); H 1 ), weak-* in L°°((0, T); L 2 ), 
and strongly in L 2 ((0, T); L 2 ). 
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It is straightforward to check (e.g. [TcmOlJ, section 3.2) that addi- 
tionaly 

Um-V<^u-Vu in L 1 (Qt) 

u m • V <i, u m " Vu 4 U • Vt>, U • Vu> in L^Qt), 

«V^, Vu) -» (uV^, Vi/), «V^, Vi/) -> (t«V0, Vi/), 
inL^O,? 1 ) Vi/ e H 1 . 

One easily verifies that (u, v, w, <f>) is a global weak solution of (fT|)- (fTT|) 
(c.f. [TemOlj . Chapter 3.) 

Let f = Acj)V(j) and = A^V^. By Fubini's theorem and the 
above, the functions v e m ^w e m converge to v,w, in the strong topology 
of L 2 (Qy). By regularity of solutions to the Poisson equation, V<^ <g> 
V<fim converges to (g) V</> in the strong topology of L 2 (Qt). Using 
the characterization (|39p and the estimate (|60p . one verifies through 
multiplying by test functions that 



u is a weak solution of (|22p . Let < s. By proposition ([2]) we 
conclude u t G L 2 ((s,T);H 1 ) n L°°((s, T); L 2 ). Applying the stan- 
dard regularity results for the Navier-Stokes system to (|22p . we find 
u G L°°((s, T); H 2 ). We conclude from the Sobolev embedding, 



for some a > 0. By choosing p = 2, proposition [T] implies that v,w £ 
L 00 (Q Sj r). The C 2+a regularity u and w follows as in proposition 



f GL°°((0,T);V*), i t GL 2 ((0,T);V*). 




u t | b + |Vu| b dxdt < oo. 



This is sufficient to conclude that 



u G C Q (Q S , T ) 
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[H although here the boundary data on the bottom of Qt may not 
be smooth. By [Lie96j , Theorem 6.44 v,w G C^Q^t)- Then V0 G 
C a (Q Sj r) as well and from the material in [Lic96], Chapter 5, v,w G 
C 2+a (Q Sj T)- Serrin's result (c.f. |Ser62j). which states that a weak 
solution of the Navier-Stokes equations in two dimensions with Ut G 
L°°((s, T); L 2 ) and f G C a (which is indeed the case here), enjoys 
u G C 2+Q (^). 

5 Entropy and Decay 

In |BD00j . Biler and Dolbeault studied various convergence estimates 
for weak solutions of the Debye-Hiickel system to the static solution. 
The trick was to observe that in certain instances entropy produc- 
tion terms can be bounded in terms of a relative entropy. The same 
observations can be made in the case of the system modeling electro- 
hydrodynamics. Theorem [3l which is modeled after the work of Biler 
and Dolbeault, follows by proving that the relative entropy of the 
electro-hydrodynamic system decays atleast with a rate depending 
only on £1. 

Let v,w G L 2 and <fi G Hg be a solution of the Poisson equation 
©. Let (voo, Woo, 0oo) be the stationary solution from Corollary [TJ 
Define the relative entropy 

W" re i = y ^oc(w) +4>w x (w)dx + -||V0- V^oo|| 2 2 + -||u|| 2 2 . 

Assume is bounded and uniformly convex subset of R n . We will 
show that there is a constant Ai depending only on Q so that that if 
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(u,v,w,<f>) is global weak solution of (fTj)- (fTT|) . then 

^<-XW veh Vt 6(0,(50). (61) 

Lemma 11 (Generalized Csiszar-Kullback Inequality, [UAMTOO]). 
Let »,wgL log L and <f> G Hj. T/ien 

||« - WooIIli + W W -Woo\\ L i + ||V0- V0oo||^2 + |MIl2 < 4Ty re i. (62) 

Combining ([6"T]) with ([62]) and setting ei = 4W re i(0), theorem ([3]) 
is now proved. 

We now prove (]6ip . we need the following lemmas. 

Lemma 12. Let v, w £ L 2 and 6 Hg 6e a solution of the Poisson 
equation ([5]). Ze£ 

1 



Woo = / 4>(Voo) +lp(w OQ )dx + ^||V^>oo|Il2 



n 



2' 



and 



e v e 
vri = I v dx f . ^ - , WM = / -wodx- 



Let J be defined by ([32]) and W defined by ([T3]) . T/ien 

^ rel = W + Woo (63) 

^v M (v) + i M (i») + ^|n| 2 dx + J[^oo] - ■/[</>]. (64) 
n 1 

Proof. Equations (j63H and (|64f) follow from elementary manipulations 

involving the definition of the Maxwellians vm and wm and relations 

©, (HO} and 

v dx = / t>oo dx = / vq dx, / w dx = / dx = / wo dx. 
<> Jf2 Jo, Jn 

□ 
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Lemma 13 ([BDOQ]). Let Q, be a bounded, uniformly convex subset 
of R n . There exists C13 depeding only on fl so that for all v,w G 
LlogL, <fi & Hj, 



+^w M (w)dx 

< C13 / v|Vlog(we~^)| 2 + w|Vlog(>e^)| 2 dx. 

Proof of d61|). Note that Woo is independent of t. Thus, if (v, w, <f>, u) 
is a global weak solution of (fTj)- ([TTj) . we have, according to (fT8|) and 

HMD, 

~~ IF 

= - / u|Vlog(i;e-*)| 2 +tB|Vlog(tBe*)| 2 + |V«| 2 ds. 
Jn 

By the lemma (|13p . 

dW ^ <r I ^ ( ^ ( W \ A IIV7 ||2 
< -C13 / ^log + wlog dx - \\Vu L 2. 

at J Q \vmJ \wmJ 

By the Poincare inequality, there is C depending only on Q for which 

—7— < -C13 / flog + wlog dx - — \\u L2 . 

fit* V^M/ \W M J 2 

Let Ai = min{ci3,C}. Then, applying (fM|) . 

^<-Ai(W rcl + JM- J[<U). 

Since is a minimum of J (see proof of theorem [1]), the difference 
of the last two terms is positive. This gives (|6ip , □ 



6 The 3 Dimensional Case: Small Data 

For domains with more general geometry, the techniques from the pre- 
vious section are difficult to apply. The main reason for the difficulty 
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is that W re \ is not quadratic and it is not clear to what extent the 
logarthmic Sobolev inequality of lemma (|13p depend in the domain 
geometry. 

To remedy this difficulty, we study a linearization of the relative 
entropy W Te \ about the stationary solution (Voo, Woo, (poo)- We are able 
to show that the linearization satisfies a decay estimate similar to 
(fUTj) . This approach is succesful because W ve \ is locally quadratic 

about (Uoo.Woo^oo). 

In order to construct the linearization of the relative entropy W re \ 
about the stationary solution {voo, Woo, (poo), consider an expansion of 
ip r (s) = s log(s/r) — s + r for s > 0, r > 0. By Taylor's theorem, 

A(s) = tp r (r) + (s - r)Vv(r) + -(s - rfip"(r) + 0((s - r) 3 ) 

If s = v(x,t) (resp. w(x,t)) and r = Voo(x) (resp. Woo(x)), the leading 
order term in v,Voo,w, Woo of the integrand of W Te \ is 

1 (V ~ Voo) 2 l (w- Wqq) 2 

2 Voo 2 Woo 

Motivated by this observation, we define 

Jn 2 2voo 2voo 

(65) 

Lemma 14 (Weighted Poincare Inequality). Let Q be a connected, 
open subset of W 1 and < a < b < oo. Then there exists C14 = 
cn(a,b,Cl) such that if p G H 1 satisfies 

o < p(x) < b, a.e. x £ Q 
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and f £ H 1 satisfies 

fdx = 0, 



J 



m 

then 



[ f 2 dx<c 1A [ \V(fp)\ 2 dx. 
Jn Jn 

If p is merely positive and p^ 1 is integrable, then there is cu = 
cu(p, O) for which the same conclusion holds. 

Proof. Suppose that no such constant exists. Then there is a sequence 
of functions {fi}^ and {p;}^ in H 1 ^) with 

/ f t dx = 0, \\fi\\l2 {n) >i [ \V(f tPl )\ 2 dx. 
Jn Jn 

Let hi = /i/||/j|| L 2(Q) so that 

iHIMi^n) >* / \V{hiPi)\ 2 dx. (66) 



Write <7j = hiPi. Then 



n Jn 



gf dx= h^ p^ dx < b 2 



shows that gi is bounded in H 1 (0) and thus converges weakly to an 
element g € H 1 (Q). Fatou's lemma, 

/ Wg\ 2 dx < liminf / \Vgi\ 2 dx< lim - = 
Jn Jn i 

and the connectedness of Q shows that g(x) = G a.e. for some constant 

G. Moreover, p^ 1 are uniformly bounded in L 2 and so we may extract 

a subsequence (reindexed by i) converging weakly to a in L 2 with 

b^ 1 < a < a -1 a.e. Clearly, 



lim / p i l gi dx = I ag dx. 

*->°°./n Jn 
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Then 

G a dx = / ag dx = lim / p~^ gidx 

= lim f hidx = lim --— f fidx = 0. 

i-^VJCl ^° ll/i|lL 2 (Q) Vn 

We infer G = since is nonzero on a set of positive measure. The 
contradiction with gives the existence of C14 = 014(0,6,0). 

The existence of C14 = c\±(p,£l) follows from the above proof by 
setting Pi = p for each i. □ 

Lemma 15. Let 0, be a bounded subset o/R n with smooth boundary. 
Then there exists positive constants C15 , c\q and cyj depending only on 
Q, Voo and Woo such that if (u, v, w, <f>) is a weak solution of (fT1)- (fT0l) 
on Qt then 

^7 < -C15L + Ci 6 L 2 + Cl 7 [| V(0 - 0oo)Hl2, 

dt (67) 

a.e. t G (0,T). 

Remark 5. The following proof applies equally well to the modified 
Galerkin approximation from lemma [31 

Proof. For the sake of clarity, we will assume that w = Woo = 0. The 
general case requires only minor modifications. Define 

e = v-w o0 , ^ = ( ^_0 oo . 

Using ([5]) and ([30]) . we compute 

dL f e „, „, , 
— = / u u t H e 4 + • Vipt dx 

= / -|Vu| 2 + u • VcfiAcfidx + / (— if) ) v t dx = A + B. 

Jn Jn \Voc ) 
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Write 

Ai= -|Vu| 2 dx, A 2 = u-X7(pA<pdx. 
Jn Jn 

Note that by corollary Q3 Vfoo = fooV^oo. Note also that 



Vv = Voo V ( — ) + wV^oo. (68) 



Using the fact that v is a weak solution, 

B = - / V ( — - ip ) ■ (Vv - vV6 - vu) dx. 

Jn \Voc J 

Then, computing with ([68]) . 



B 







/ ^oo 




In 


\VooJ 



2 



dx 



eV ( — ) • Vtp + 2^ V ( — ) Nipdx 

n \VooJ \Voo. 



v\Vi)\ 2 dx 

n 



+ J UN I — J - uVV'j • udx 
= Si + B2 + -B3 + S4. 

We proceed by bounding ^2, -62,-83 and -B4 in terms of A\ and i?i 
and integrals of higher powers of |e| and |V?/>|. 

Remark 6. The presence of transport of the charges by the velocity 
make the following calculation somewhat more subtle than the analo- 
gous analysis for the Debye Huckel system, c.f. [BAMV04]. As will be 
demonstrated immediately below, the net exchange of kinetic energy 
tj|u| 2 and the relative energy ^f— is a second order contribution. This 
is in agreement with the cancelation of the entropy production due to 
transport and the kinetic energy production due to forcing seen in the 
derivation of the basic energy law (|18p 
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We have 

A 2 + B A = J u • ^A0V(/> + A0V (j-^j ~ A0V^ da. 
Adding and subtracting A^ooV (^r~l an( ^ using the relation v 

A0, Voo = A^oo, 

A 2 + B 4 = / u- ( A^V^oo + A0ooV ( — ) +eV( — J J da;. 
Finally, using the relation (|68|) . 



vl 2 + 5 4 = y u |^Vw + eV I — J j dx. 

Because u is divergence free and so is orthogonal to Vi> in L 2 , the first 
product in the integrand vanishes. By Young's inequality, 

A 2 + B A < -aiBi + a 2 / |u| 2 e 2 dx (69) 

Jn 

where a\ and a 2 are positive and a\a 2 = 4. 
The term B3 is nonpositive. 

To estimate B 2 we first note that Voo is bounded. Thus there exists 
a C 2 depending only on for which 

B 2 < -(h + 2d 1 )B l + b 2 [ e 2 \ViP\ 2 dx + 2C 2 d 2 [ |VVf dx. (70) 

Jn Jn 

where b\, b 2 , d\,d 2 are positive and b\b 2 = d\d 2 = 4. 

Let p = — and / = e. By corollary [IJ p satisfies the first hypoth- 
esis of lemma HH and the integral of e is zero, satisfying the second 
hypothesis of the lemma. Hence there is a constant C14 depending 
only on and f2 for which 



e 2 da; < C14 



n 



2 

dx. 
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/ 

Jn 



2e 2 



3 is C 3 








/ Voo 


H-) 


Jn 


\VooJ 



2 

dx = -C 3 Bi. (71) 



In v o 

Similarly, since tp G Hq is a solution of the Poisson equation with 
right hand side e, there is also a C4 = 6*4(^00, O, C14) for which 

/ ^|VVfdx< -CaBl (72) 
2 

Finally, by the Poincare inequality, there is C5 = Cs(f2) for which 

f h u \ 2 dx < -C 5 Ax. (73) 
Jn 1 

Adding ([7"l~|) .(fT2 ]) and ([73]) together, we have shown that there is a 
positive constant Cq = Cq(voo, f2, C14) for which 

L < -CfeOAi + Si). (74) 

Arguing in a similar fashion, there is clearly a C7 = 6*7(0, i^, a 2 , 62) 
for which 

/ a 2 \u\ 2 e 2 + b 2 e 2 \V^\ 2 dx < C 7 L 2 . (75) 

By (|M]),{IZnD and JSJ, we have 

^ = Ai + A 2 + #i + #2 + B 3 + £ 4 

< Ai + (1 - a x - h - 2d 1 )B 1 + 2C 2 d 2 [ | V?/f dx + C 7 L 2 . 

Jn 

Choose 01, 61 and di so that ai—bi—2d 2 = |, thereby fixing 02, 6 2 , <i 2 , C7. 
Then, using (JUD and the fact that Ai + Bi< 0, 

1 



n < --^L + C 7 L 2 + 2C 2 d 2 [ \V^\ 2 dx. 
dt 2Cq 



Setting C15 = 2<Iv Cl6 = ^7 an< ^ Cl7 = ^C 2 d 2 , the lemma is now 
proved. □ 
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From theorem Q] there is p\ > so that for all < po < Pi, 

|0oo(^)| < 1, Va€fi. 

Following the proof of lemma (|15p and using corollary [H one may keep 
track of the constants Ci5,c%e and cyj to prove 

Corollary 2. There exist 'positive constants c' 15 , d ie , c' 17 and pi de- 
pending only on ft such that if po < p\ then 

— < -c' 15 L + p c' 16 L 2 + p c' 17 \\X7{(f> - </>oo)Hl2, 

V* € [0,T ]. 

6.1 Proof of theorem [4] 

We are assuming Q C R n , n = 2, 3 is bounded with smooth boundary. 
We will first use lemma [T5l to prove an extension property analogous 
to theorem [6) 

Let uo € H,uo,iOo € L 2 . Let {vq,Wq} be a sequence of functions 
satisfying with , it>o "^o in L 2 as /i j 0. Let (u^, u^, w^, </>^} 

be the local modified-Galerkin approximate solution on Qt obtained 
from lemma [3] with initial data uo , Vq , Wq . 

Let p\ be the constant from corollary [2l Since ||V($J, — </>oo)Hl2 < 
2L, we may choose 

so that if po < m in{pi,p2} then 

— < -^L + p c' 16 L 2 , a.e. i e [0,T ]. 
(see remark O) Finally, let 

4 P0Ci 6 
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Then 



L(0) < 5 implies L(t) < be 



a.e. t G [0, T ] 



(76) 



This inequality implies that the L 2 norm of u^, and w m remains 
bounded independently of t. The extension property now follows ex- 
actly as in theorem [U One checks (just as in the proofs of estimates 
([59]) and (f60|) ) that there are constants ci8,cig independent of m and 
h for which 

l-T 



SUp ||u m , D m , U7 m ||L2 

te(o,T) 



J o 



and 



where 



+ 



sup 

te(o,T) 



dt 



| f /ii|2 , 



+ 



c?t dt 



ds < c\s- 



H 1 



dill 



dt 



ds < Cig 



Letting — > 0, m — > oo we see that some subsequence of (u 



h n ,h 



h ih 



IT 111 . 

mi ^mi " J m'> "ml 

converges to a global weak solution (u, v, w, (ft) of (fT])- (fTTj) . This proves 
the first part of theorem HI 

Note that S 2 = 2L. From £[6]), let e 2 = 25. If <f 2 (0) < e 2 , then 
L(0) < 5 and we may apply (f76|) for a.e. i G [0, T]. Let A 2 = ^ which 
depends only on $7. Then the global weak solution (u, v, w, 4>) satisfies 



&> = 2L < 25e 



-t\ 2 



a.e. t G [0,T]. 



The last part of the theorem now follows immediately from the 
following lemma. Using the embedding V C H, minor modifications 
of the proof of [BHN94J , theorem 6 gives 
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Lemma 16. Let £1 C M n ,n = 2,3 be a bounded, open set with smooth 
boundary. If (u, v, w, (j>) satisfy ([5]), (fTUj) . 

SUp \\v, ||l 2 < DO) 

te[o,oo) 

sup / -ulVlog^e"^)! 2 + wlVlog^e*)! 2 + |Vu| 2 dx < oo, 

te[0,oo)JQ t 

then there is some sequence tj — ► oo for which 

lim L(tj) = 0. 

6.2 Proof of theorem 

Let u G V n H 2 , vq, wo £ H 2 and T > 0. Let {vq , } be a sequence 
of functions satisfying (pl2"j) with Vq,Wq — ► vq,wo in H 2 as fa j. 0. Let 
(u^j, u^j, w^, 4>m) be the local modified-Galerkin approximate solution 
on Qy obtained from lemma [3] with initial data uo, v$ , U)q . Let p2 and 
e2 be the constants from theorem HJ Assume that po < p^ and <^ < £3 
where p% < P2 and €3 < £2 will be determined below. Applying the 
results from theorem [H (u^ n ,v!^ n ,w^ n ,4>^ n } is defined for all t G [0, T] 
and some subsequence converges to a global weak solution of ([!])- 
(fTT|) as m — > 00 and fa — ► 0. In the sequel we suppress the sub- and 
superscripts m and h. 

From lemma [3] we infer that vt and wt are smooth in x. In partic- 
ular, vt and tot are classical solutions of 

v tt + u t ■ Vv + u • Vv t = V • (Vu t - v t V<p ~ vV<p t ), 
w u + ut ■ Vw + u • Vwt = V • (Vw t + w t V(j) + wV4> t ). 

Multipling these equations by vt and wt and integrate over O. Inte- 
grating by parts and noting that the boundary terms vanish in this 
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case as well, 



1 d 

2dt 



Oi 



Of 



+ 



L 2 



„dv 



where 



du „dd> „ , dv 



We estimate the L 2 norm of <£. From 0, the integral of v$ vanishes. 
By the Sobolev embedding H 1 C L 4 and the Poincare inequality, the 
estimate 



l$l 



L 2 



< c 



„9u 
V — , V — 

at' at 



+ 



2 

L 2 , 



for some C = C(f2) is straightforward. The analogous estimate holds 
for 

Now differentiate (|47p with respect to t and multiply the resulting 
system componentwise by U{ and add the equations for % = 1, . . . , m. 
One concludes 



2 eft 



u 



at 



+ 



L 2 



„<9u 



2 

L 2 



du du 



dt' 1 dt 



b[—,u,— ) + *,V— 



7 dv\ 



dt J 



where 



dt dt 



The usual estimate in the small data regularity proof for Navier-Stokes 
shows that 



, 'du du\ 



dt 



iVul 



L 2 - 



Using the embedding H 1 <^-» H 1 from the Poisson equation, there is 
C = C(J2) for which 



1^1 



L 2 



< c 



„dv „dv 
V — , V — 
dt' dt 



\v,w\\ ul 
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Define 



G(t) = 

H(t) = 
I(t) 



du 


Ov 


dw 


2 




dt 


~dt 


L 2 


du 


dv 


dw 


2 




w 


~dt 


H 1 



\u,v,w\\ u i. 



The above estimates and the Poincare inequality show that 
d 



dt 



G(t) + H(t)(l-I(t))<0, VtG[0,T]. 



Now we develope a relationship between I(t) and G(t). From ([47|) . 

du 



2||Vu[|£a = -2(V0 8) V0, Vu) - 2 

< nvui|2 2 + wvm, + 



at' 



u 



du 



dt 



I u IIl2- 



Define 



Then 



C 20 = SUp ||V0||l4, C21 = SUp ||u|| L 2. 

te(o,T) te(o,T) 



||Vu|| H i < c 20 + c 2 iG(t)2, Vie[0,T]. 
Similarly, by (JHj) and the triangular inequality, 



MIhi < 2 ( / v dx ) + 3||V^||£ 2 



2 Po-3( ^,«) +3(vV^V« 



< 2pg + 3 



9* 
(9t> 



5/ 



L- 



H| L 2 + CiC 2 o|MIhi 



for some Ci = Ci(fi). Applying a similar estimate the it?, we find 

1 



1 - C1C20 



4/9 + C22 
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cJu dw 
~dt , ~dt 



where 

c 2 2 = 3 sup \\v, w\\ui. 
te(o,T) 

Define 

4pp , c 22 

C23 = 1 ~ hC 20 , C 2 4 = 7^ hC21. 

1 - C1C20 1 - C1C20 

We have shown that 

1(f) < c 23 + c 24 G(t)i 

Thus, 

^G(i) + ff(t) (l - c 23 - c 24 G(t)l) < 0, Vt e [0, T]. 

Assume for the moment that p 3 , 63 and 63 can be chosen so that 

l-c 2 3-C2 4 G(0)5 >0. (77) 

independently of h and m. If follows that 

G(t)+ f H{s)ds<G{0), te[0,T] (78) 
Jo 

for all m and h. Letting h — ► and m — > 00 we may extract a subse- 
quence of (u^j, u^j, u;^, which converges to a global weak solution 
(u,v,w,4>) of (fB- (fTT1) . From (f78l) . this solution satisfies the estimate 

uGL°°((0,r);H 2 ), u t GL°°((0,r);L 2 ). 

Arguing as in the end of the proof of theorem ([2]), we find that u is 
Holder continuous on Q T and C 2+a on compact subsets of Qt and 
v, w are C 2+Q on Q s ,t for any < s < t < T. The standard arguments 
show that (u,v,w,(/>) is unique, c.f. (TemOlj . 

Now we show that /J3, £3 and 83 may be chosen in order that (|77p 
be satisfied. Note that ([77]) holds provided po, 020,021 and C22 are 
sufficiently small and G(0) is sufficiently small with respect to 22a. 
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By the regularity of solutions to the Poisson equation, there is 
C 2 = Cfe(n) so that 

C20 < C2c| 2 - 

By the triangular inequality, 

C22 < sup 3{||v - «oo,^- ^oo||l2 + ||«oo,^oo||l2} 
te[o,T] 

< sup C 3P0 {VL + 1) 

te[o,T] 

for some C3 = C^(U), provided po < p\. Finally, 

C21 < 2 sup VX. 

te[o,T] 

By assumption, po < P3 < P2 and <f2 < £3 < £2- By theorem HI we 
have then that 

sup L = sup -Si < e 3 . 

te[o,T] te[o,T] 2 

Using the above bounds on C2o,C2i and C22 in terms of L and the 
bound on L in terms of £3, we may choose €3 and p% so that 

C23 * i-dc^pgfa + i) + C2Clpl ^ + 1)<1 r 

g 3 po(^i+l) 
C24 -l-6C 1 C 2 C 3 po(^+l) +2 ^ <1 ' 

G(0) < \ 

This implies (|77j) . Certainly there exists C4 = C^{Q) so that 

G(0) < C 4 ||uo, uo,^o|Ih 2 - 
Setting J3 = Ijt- now completes the proof. 

V 4C/4 
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7 Conclusion 



The equations of a viscous, incompressible fluid coupled with diffuse 
charges in two and three dimensions have been studied. The key step 
toward the existence of global in time solutions is the presence of a 
decaying entropy function which guarantees the dissipation of kinetic 
and electrostatic energy and entropy. 

The most serious obstruction to formulating a global existence re- 
sult like theorem [2] when dimfi = 3 is the Debye-Hiickel system. A 
different approach is to consider very weak solutions, i.e. those where 
v,w £ L°°((0, T); L 1 ) satisfy the weak formulation in terms of test 
functions u £ C 1 ^). How one defines the forcing term Ac/>V0 and a 
solution of ([5]) then becomes a more delicate matter. 

The techniques used in this paper certainly apply to other Dirich- 
let conditions than (llOj) and other 2nd order elliptic operators than 
the Laplacian. A future avenue of study are electrorheological fluids 
where the charge is vector valued and the potential is the polarization 
potential, see [ZGL + 08] . 

References 

[AMTU01] Anton Arnold, Peter Markowich, Giuseppe Toscani, and 
Andreas Unterreiter. On convex Sobolev inequalities and 
the rate of convergence to equilibrium for Fokker-Planck 
type equations. Comm. Partial Differential Equations, 
26(l-2):43-100, 2001. 



58 



[BAMV04] Naoufel Ben Abdallah, Florian Mehats, and Nicolas 
Vauchelet. A note on the long time behavior for the drift- 
diffusion-Poisson system. C. R. Math. Acad. Sci. Paris, 
339(10):683-688, 2004. 

[BD00] P. Biler and J. Dolbeault. Long time behavior of solutions 
to nernst-plank and debye-hiickel drift-diffusion systems. 
Ann. Henri Poincare, 1:461-472, 2000. 

[BHN94] P. Biler, W. Hebisch, and T. Nadzieja. The Debye system: 
existence and large time behavior of solutions. Nonlinear 
Analysis T.M.A., 23:683-688, 1994. 

[DH23] P. Debye and E. Zurek Hiickel. Zur theorie der electrolyte, 
ii. Phys. Zft, 24:305-325, 1923. 

[EG92] Lawrence C. Evans and Ronald F. Gariepy. Measure the- 
ory and fine properties of functions. Studies in Advanced 
Mathematics. CRC Press, Boca Raton, FL, 1992. 

[FG09] Jishan Fan and Hongjun Gao. Uniqueness of weak solu- 
tions to a non-linear hyperbolic system in electrohydrody- 
namics. Nonlinear Anal, 70(6):2382-2386, 2009. 

[GT01] David Gilbarg and Neil S. Trudinger. Elliptic partial dif- 
ferential equations of second order. Classics in Mathemat- 
ics. Springer- Verlag, Berlin, 2001. Reprint of the 1998 
edition. 

[Li09] Fucai Li. Quasineutral limit of the electro-diffusion model 

arising in electrohydrodynamics. J. Differential Equa- 
tions, 246(9):3620-3641, 2009. 



59 



[Lie96] Gary M. Lieberman. Second order parabolic differential 
equations. World Scientific Publishing Co. Inc., River 
Edge, NJ, 1996. 

[RLW06] R. Ryham, C. Liu, and Z.-Q. Wang. On electro-kinetic 
fluids: one dimensional configurations. Discrete Contin. 
Dyn. Syst. Ser. B, 6(2):357-371, 2006. 

[RLZ07] R. Ryham, C. Liu, and L. Zikatanov. Mathematical mod- 
els for the deformation of electrolyte droplets. Discrete 
Contin. Dyn. Syst. Ser. B, 8(3):649-661, 2007. 

[Ser62] James Serrin. On the interior regularity of weak solu- 
tions of the Navier-Stokes equations. Arch. Rational Mech. 
Anal, 9:187-195, 1962. 

[Tay97] Michael E. Taylor. Partial differential equations. Ill, vol- 
ume 117 of Applied Mathematical Sciences. Springer- 
Verlag, New York, 1997. Nonlinear equations, Corrected 
reprint of the 1996 original. 

[TemOl] Roger Temam. Navier-Stokes equations. AMS Chelsea 
Publishing, Providence, RI, 2001. Theory and numerical 
analysis, Reprint of the 1984 edition. 

[UAMT00] Andreas Unterreiter, Anton Arnold, Peter Markowich, 
and Giuseppe Toscani. On generalized Csiszar-Kullback 
inequalities. Monatsh. Math., 131(3):235-253, 2000. 

[ZGL+08] J. Zhang, X. Gong, C. Liu, W. Wen, and P. Sheng. Elec- 
trorheological fluid dynamics. Physical Review Letters, 
101(194503), 2008. 



60 



